Strain Relief in Cu-Pd Heteroepitaxy 
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We present experimental and theoretical studies of Pd/Cu(100) and Cu/Pd(100) heterostruc- 
tures in order to explore their structure and misfit strain relaxation. Ultrathin Pd and Cu films are 
grown by pulsed laser deposition at room temperature. For Pd/Cu, compressive strain is released 
by networks of misfit dislocations running in the [100] and [010] directions, which appear after 
a few monolayers already. In striking contrast, for Cu/Pd the tensile overlayer remains coherent 
up to about 9 ML, after which multilayer growth occurs. The strong asymmetry between tensile 
and compressive cases is in contradiction with continuum elasticity theory, and is also evident in 
the structural parameters of the strained films. Molecular Dynamics calculations based on classi- 
cal many-body potentials confirm the pronounced tensile-compressive asymmetry and are in good 
agreement with the experimental data. 

PACS numbers: 68.55.-a, 68.35.Gy, 81.15.Fg 



X3 
O 

o 



> 
oo 

m 
m 
o 

o 



o 



X 
S3 



Heteroepitaxy with different thicknesses ranging from 
monoatomic layers up to micrometers produces various 
structures, and is one of most important routes to artifi- 
cially obtain functional materials and devices. Stress and 
its relaxation due to an interface lattice mismatch are an 
essential problem in the heterostructures. Classical con- 
tinuum theory defines the concept of an equilibrium crit- 
ical thickness above which the lattice stress relaxes 
with the introduction of a misfit dislocation (MD) 
The critical thickness is determined by energy balance 
between strain energy buildup and strain relief due to 
dislocation nucleation in the mismatched structure. The 
value of he predicted by the continuum theory, however, 
often does not agree with experiments. Moreover, the 
critical thickness in the continuum elasticity theory is in- 
dependent on the strain type, be it tensile or compressive. 
Atomic details of the interface structure such as surface 
steps and surface roughness, are usually ignored. Actual 
stress relaxation of coherent strained films is kinetically 
limited, typically with high barriers. Recently, atomistic 
studies on model systems have been done to determine 
the transition paths and energy barriers from coherent to 
incoherent states in real epitaxial films ,^^, .4]. One of 
the main findings has been a striking asymmetry between 
the tensile and compressive cases, due to the anharmonic- 
ity and asymmetry of the atomic interactions 0, 0| . Fur- 
ther, the concept of a size-dependent mesoscopic mis- 
match has been proposed to explain strain relaxations in 
the early stage of homo and heteroepitaxial metal growth 
on the atomic scales 

In this Letter, we present results from a combined ex- 
perimental and theoretical study on a model system of 



Cu-Pd heteroepitaxy grown by pulsed laser deposition 
(PLD). The bulk lattice parameters of Cu and Pd are 
flcu — 3.6lA and a-p^ = 3.89A, respectively. The lattice 
misfit induces a large compressive (to = —7.2%) strain 
in the Pd overlayer on Cu(lOO) and a tensile (m — 7.8%) 
strain for the Cu overlayer on Pd(lOO). The Pd films 
on Cu(lOO) and the Cu films on Pd(lOO) were grown 
by pulsed laser deposition (PLD) in a multi-chamber 
ultrahigh vacuum (UHV) system with a base pressure 
p < 5x 10~^^ mbar, andp < 2x 10~^° mbar during depo- 
sition. Prior to deposition the copper and palladium sub- 
strates were cleaned by cycles of Ar+ sputtering followed 
by annealing at 873 K for Cu substrate or 950 K for Pd 
substrate until clean Auger electron spectroscopy (AES) 
spectra, sharp low-energy electron diffraction (LEED) 
spots, and atomically smooth terraces under scanning 
tunneling microscopy (STM) were observed. The sub- 
strate temperature was kept at room temperature during 
deposition. The substrate was placed about 100 ~ 130 
mm away from targets. A KrF excimer laser beam (wave- 
length 248 nm, pulse duration 34 ns, typical pulse energy 
300 — 350 mJ and pulse repetition rate of 3 — 5 Hz) was 
focused onto targets. During deposition the growth pro- 
cess was monitored by a reflection high energy electron 
diffraction (RHEED). All STM measurements were per- 
formed in the constant current mode at a 0.2 — 0.5 V 
positive tip bias and a 0.1 — 0.5 nA tunneling current. 
Compared to thermal deposition (TD), an extremely high 
instantaneous flux of atoms in PLD favors much larger 
nucleation density 0, Q • This helps to suppress the for- 
mation of Cu(001)c(2 X 2)-Pd ordered alloy in the initial 
growth stage for Pd/Cu(100) ^ and the phase transition 




Figure 1: STM images for (a) 0.5 ML Pd (200 x 200 nm^); (b) 
1.6 ML Pd (200 X 200nm^); (c) 4.8 ML Pd (100 x 100 nm^); 
(d) 6.0 ML Pd (100 x 100 nm^ deposited at 300 K. 

at 1 ML from fee to bet strueture for Cu/Pd(100) 
Even though the present STM setup does not provide a 
direct distinction of the Pd and Cu atoms, neither pits 
in the substrate surface nor eroded substrate step edges 
show any evidence of Cu-Pd interlayer mixing. Moreover, 
the surface area covered by the islands for submonolayer 
coverages, as determined from STM images, are always 
consistent with the coverage from RHEED. The RHEED 
data confirm the layer-by-layer growth mode for both 
systems. It persists up to at least 5 — 6 ML in Cu-Pd 
heteroepitaxy. 

Figures ^a)-(d) display STM images taken at room 
temperature for Pd films on Cu(lOO). For submono- 
layer Pd at 0.5 ML, many islands with a monolayer 
height of about 1.80 A were observed as shown in the 
inset of Fig. ^a). With increasing Pd thickness such 
islands are always formed on a completely filled under- 
layer up to a total thickness of about 3 — 4 ML, indicat- 
ing ideal 2D growth mode. Above 4 ML a striking fea- 
ture is the appearance of misfit dislocations as marked 
by many protruding stripes with an average height of 
Ah = (0.57 ± 0.10)1 (see Fig. [He) and (d)), by which 
the compressive stress in Pd film relaxes efficiently. They 
align nearly regularly along the [100] and [010] directions. 
The average spacing of MD lines is (7.09 ± 0.30) nm for 
4.8 ML. The layer-by-layer growth persists up to 6 ML, 
after which multilayer growth mode occurs. 

The corresponding layer-by-layer growth mode for 
PLD Cu films on Pd(lOO) has also been confirmed by 
STM images as indicated in Fig. |21 At 0.85 ML, sec- 
ond layer nucleation is not found, while in the 2.1 ML 
thick film more than 95% of the second layer is accom- 
plished (Figs. My>) and (c)). At 2.1 ML, few of the is- 
lands in the third layer have a preferential orientation. 



whereas at 3.05 ML most of the fourth layer islands are 
rectangular along the (110) directions. With increasing 
Cu thickness to 6.0 ML, 90% of the sixth layer is com- 
pleted, where nearly all the pits in the sixth layer and 
the seventh layer islands are rectangular. Above 6 ML 
the morphology of Cu films is changed because of the 
multilayer growth mechanism. A typical image for 9.0 
ML is given in Fig. [5Jf). Four layers appear simultane- 
ously on the surface. Moreover, only 70% of the ninth 
layer is completed while the remaining 30% appears as 
the tenth and eleventh layers. It should be pointed out 
that the height or depth of all (ir)rcgular islands and 
pits is about 1.65 A until 9.0 ML as shown in the in- 
set of Fig. [2Id) . This vertical magnitude corresponds to 
the interlayer spacing of highly strained fee Cu [l3 . In 
PLD Cu films on Pd(lOO) we have not found any trace 
of dislocation nucleation up to 9 ML. Above 9 ML, the 
multilayer island structure with the (110) orientation is 
still visible and the film surface becomes more and more 
rough. 

IV-LEED measurements of the specularly scattered 
electron beam were done to obtain structural informa- 
tion on the films. The average interlayer spacing of 
PLD films was calculated from LEED intensity of the 
(00) diffraction beam vs. energy curves based on a ki- 
netic model Fig. Ogives the layer thickness depen- 
dence of the interlayer distance for both systems. In the 
case of Pd/Cu(100) films, the interlayer distance is about 
1.82(5) A up to 3 ML, very close to the value of Cu sub- 
strate, and increases rapidly above 3 ML to reach the 
value of bulk Pd. For Cu on Pd(lOO) up to 10 ML, the 
value stays constant (« 1.64 A). This agrees well with the 
measured height of monolayer islands by STM. Above 10 
ML the interlayer spacing gradually increases close to the 
value of bulk fee Cu at a coverage of 20 ML. 

The actual in-plane lattice constant of deposited films 
was obtained in situ by measuring the spacing between 
the reciprocal lattice rods in the RHEED image. The 
RHEED spacing of substrates serves as a reference. The 
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Figure 3: Perpendicular interlayer distance calculated based 
on a kinetic model from LEED intensity of the (00) diffraction 
beam vs. energy curves. The change of the interlayer distance 
occurs at about 3—4 ML for Pd/Cu(100) films and at about fO 
ML for Cu/Pd(100) films. Assuming a pseudomorphic growth 
up to 3 ML for Pd on Cu(lOO) and 10 ML for Cu on Pd(lOO), 
PLD films can be identified as fee Pd and fct Cu structures. 
The results from experimental and theoretical work at 300 K 
are indicated by open and solid symbols, respectively. The 
lines are guides to eye. 



in-plane lattice strain is defined as (ab„ik — af„„,)/abuik(%). 
Fig. 01 gives the measured data. For Pd/Cu(100), the 
measured in-plane strain exhibits a gradual decrease with 
increasing thickness of the Pd film. The compressive 
strain is nearly fully relaxed at 9 ML. In contrast, the 
in-plane tensile strain shows a very different thickness de- 
pendence in the case of PLD Cu films on Pd(lOO). Until 
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Figure 4: In-plane lattice strain derived from the separation 
of the first order rods in the RHEED pattern as a function 
of Pd or Cu coverages. Our results from experimental and 
theoretical work at 300 K are indicated by open and solid 
symbols, respectively. See text for details. 




Figure 5: Typical configuration of stacking faults for 7 ML 
Pd on Cu(lOO) based on molecular dynamics annealing. The 
lateral size is 20 x 20 atom^ . The misfit dislocations orientate 
along the (100) direction of the substrate. 



9 ML only 0.6 % reduction in the tensile strain is found. 
Above 9 ML a step-like decrease is visible and a residual 
strain is observed until 20 ML. 

The most striking discovery in our data is the strong 
asymmetry in the strain relaxation behavior between the 
tensile and compressive strained films. This result is in 
direct disagreement with the continuum elasticity theory, 
which predicts a symmetrical behavior for tensile and 
compressive systems. 

In any growth mode, there is always the question 
whether the final configuration is dictated by kinetic con- 
siderations or it corresponds to a true thermodynamic 
equilibrium state. Given the lack of precise characteri- 
zation of the PLD growth conditions, a full theoretical 
modeling of the kinetic growth process under the experi- 
mental conditions is unfeasible. Instead, we have chosen 
to perform equilibrium molecular dynamics simulations 
of the epitaxial layer after deposition. The object here is 
to study the thermal excitation and nature of dislocations 
in the epitaxial film and to see whether they can account 
for the observed stress relaxation and tensile compres- 
sive asymmetry in Cu-Pd heterostructures. Our atom- 
istic model contains five layers of substrate and varying 
number of film layers. Two bottom layers of the sub- 
strate were fixed. Periodic boundary conditions were ap- 
plied in the plane parallel to the surface. The lateral size 
of the systems studied varied between 20 x 20 to 40 x 40. 
Interatomic forces between particles in the system were 
computed using Embedded Atom Model (EAM) poten- 
tials for Cu and Pd The initial configuration of the 
film has adatoms occupying fee positions in registry with 
substrate. 

For Pd/Cu(100), after thermal equilibrium is reached 
at 300 K, the system relaxes and gains partial strain re- 
lief through generation of misfit dislocations starting at 
2 ML and above. To facilitate a direct comparison with 
the experimental data in Fig. |31 we computed from our 
simulations the average interlayer distance for the top 
layers of the film (up to three to take into account the 
surface sensitivity of LEED). We also computed the av- 
erage in-plane lattice constant for the top three surface 
layers, and then evaluated the in-plane strain using the 
bulk value of the adsorbate as a reference. The depen- 
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dence of these values on the film thickness is shown in 
Fig. ^ There is a good qualitative agreement between 
the theoretical results with the experimental observation 
showing continuous relief of in-plane strain through the 
misfit dislocations. The exact nature of the misfit dislo- 
cations is rather complicated, varying both as a function 
of film thickness and lateral size. In Fig. [SJ we show a 
typical configuration for a film of 7 ML thickness and 
lateral size 20 x 20 obtained by molecular dynamics an- 
nealing through heating to 500 K, followed by cooling to 
K. This shows a Pd film with misfit dislocations visi- 
bly aligned along the (100) direction of the substrate, in 
agreement with the experimental observations. 

Simulation studies were also performed for 
Cu/Pd(100). Here, in agreement with the experimental 
observation, the theory shows a striking difference from 
the Pd/Cu(10G) system. The pseudomorphic strained 
fct structure for the Cu film remains stable up to 
about 9 ML. Then strain is gradually released through 
localized surface defect such as vacancies and adatoms, 
rather than the misfit dislocation mechanism observed 
in Pd/Cu(100). The theoretical results for the average 
interlayer distance of the film and the in-plane strain for 
the top surface layers as a function of the film thickness 
are shown in Figs. 01 and 0] respectively. The theory 
again is in good agreement with the experimental data, 
indicating a stable highly strained epitaxial state below 
9 ML. 

In conclusion, we have studied the stress relaxation 
mechanism in Cu-Pd heteroepitaxy both through PLD 
experiments and numerical simulations. The experimen- 
tal data show that for Pd/Cu(100), the in-plane stress 
is relaxed almost immediately above 1 ML, whereas 
Cu/Pd(100) remains coherent up to about 9 ML. While 
these observations contradict predictions from the con- 
tinuum theory, they are confirmed by our atomistic sim- 
ulations using EAM potentials for the Pd and Cu inter- 
actions. The good agreement between the experimental 
data and the equilibrium simulation study leads us to 
conclude that the observed configurations of the film at 
different thicknesses correspond to equilibrium state with 
thermally generated misfit dislocations as the strain relief 
mechanism. The microscopic nature of the misfit dislo- 
cation observed here for the Pd/Cu(100) system is rather 
complicated and does not fit into the simple model pos- 
tulated previously T3|. 

On a microscopic level, the tensile-compressive asym- 
metry originates from the asymmetry of the interatomic 
interactions. The general feature of interatomic poten- 
tials is a steeply rising repulsive part at small separa- 
tions as opposed to a very slowly decreasing attractive 
part at large separations. This strong asymmetric feature 
of the interaction is what leads to the observed macro- 
scopic tensile-compressive asymmetry. We have per- 
formed additional studies using very different potentials 



(of Lennard- Jones type) with the same general feature, 
and the results indeed remain qualitatively unchanged. 
Interestingly, we find that geometry and dimension play 
an important role in the tensile-compressive strain re- 
lief asymmetry. For the 2D model, misfit dislocations 
are favored for tensile rather than compressive strained 
systems 0, ■ Clearly, there is much more to learn in 
strain relief in heteroepitaxy through further theoretical 
and experimental investigations. 
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